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ABSTRACT 

We address the problem of line confusion in intensity mapping surveys and explore 
the possibility to mitigate line foreground contamination by progressively masking 
the brightest pixels in the observed map. We consider experiments targeting CO(l- 
0) at z = 3, Lyct at z = 7, and CII at z = 7, and use simulated intensity maps, 
which include both clustering and shot noise components of the signal and possible 
foregrounds, in order to test the efficiency of our method. We find that for CO and Lya 
it is quite possible to remove most of the foreground contribution from the maps via 
only 1% —3% pixel masking. The CII maps will be more difficult to clean, however, due 
to instrumental constraints and the high-intensity foreground contamination involved. 
While the masking procedure sacrifices much of the astrophysical information present 
in our maps, we demonstrate that useful cosmological information in the targeted lines 
can be successfully retrieved. 

Key words: cosmology: theory - cosmology: large-scale structure of universe - cos¬ 
mology: diffuse radiation 


1 INTRODUCTION 

Intensity mapping is a powerful new technique for study¬ 
ing the large-scale structure of the universe. By observing 
the large-scale fluctuations in the intensity of some chosen 
spectral line, it is possible to study a population of galaxies 
without needing the kind of dept h and sensitivity required 
to fin d each galaxy individually llMadau. Meiksin, fc Ree^ 
[TqqJ). The intensity of a line in a given volume depends 
both on the detailed astrophysical processes in the emission 
region as well as the underlying dark matter density. Thus 
intensity mapping surveys contain a wealth of cosmological 
and astrophysical information which is difficult to obtain 
using other methods. 

There are many lines that can be used for intensity 
mapping. A number of experiments either proposed or 
in progress are seeking to study the 21 cm fine structure 
fine in neutral hyd rogen, such as th e Square Kilome¬ 
ter Array (SKA) JCMill^jt^lE 20041). Low Frequency 
Array (LOFAR) (|R.ottgering et all 200fil ) . the Precision 
Array for Probing the Epoch or Reionization (PAPER) 
llParsons et aLll2 010l). the Gia nt Meterwave Radio Telescope 
(GMRT) I Paciga e t al.ll201ll) ^nd the Murchison Widefield 
Array (MWA) i Tingav et al.l l2013l) . We consider here 
three other co mmonly discussed lines: the CO (1-0) rota¬ 
tiona l li ne. llRighi. H ernande z-Monteagudo. fc SunvaevI 
I2OO8I : Ibidz et ^ I2OII : IPullen et al.l boisl : 
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drogen fine ( Pullen. Pore, fc Bock 12014 ; Gnng_et^ J 2 OM ), 
and the 157.7 /rm CII fine structure line ( Cong et al .1 120121 : 
ISilva et al.ll20lj ). Other lines which have been considered 
in the literature wh ich we do not study here in clude 
molec ular hydrogen ([C ong. Cqorav. fc SantosI l2013l) and 
Hell (IVisbah Haiman. fc Brv^ 20151 ). 

Several experiments are planned to study these lines, 
including the Carbon MonO xide Mapping A rray Pathfinder 
(COMAP) targeting CO llLi et al.l I2OI5I ). the Spectro¬ 
photometer for the History of the Universe, Epoch of 
Reionization, and Ices Explorer (SPHEREx) targeting Lya 
llDore et al.l l2014h . and the Tomographic Ion ized-Carbon 
Mapp ing Experiment (TIME) targeting CII (jCrites et al.l 
120141) . As we shall see, each of these lines has its own de¬ 
pendence on the conditions within the galaxies it is emitted 
from, and studying each line has its own set of challenges. 

One major difficulty in intensity mapping surveys is the 
problem of foregrounds. Every intensity mapping survey will 
have to deal with many types of foreground emission. Fore¬ 
grounds with continuum frequency spectra such as dust or 
synchrotron are problematic, but the removal of these fore¬ 
grounds is a well studied prob l em, es p ecially for 21 cm (s ee 
for example IPi Matteo et al.l (l2002l) : IWang et al.l ll2006l) ). 
Line foregrounds are a more difficult problem. If a spec¬ 
tral fine other than the target is redshifted into the same 
observing band, it is not easy to tell the two fines apart. 21 
cm surveys are not expected to suffer from thi s issue since 
there are so few fines at such low frequencies llCong et al.l 
However surveys in other fines will require a better 
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understanding of possible line foregrounds. As we will show 
below, CO, Lya, and CII all have potential confusing lines. 

One method which can be used to remove these fore¬ 
ground lines is to cross correlate an intensity map with an¬ 
other map in a different frequency, or with som e other tracer 
of large-scale structure (IVisbal fc Loeafeoidf . Though each 
map will have its own foregrounds, the signals from the two 
maps will be correlated, and their foregrounds will not, leav¬ 
ing behind a cross spectrum which only depends on the two 
target populations. However, this method has two issues: 
The first is simply that it requires a more complicated (and 
costly) observation, since it requires observation of two sig¬ 
nals. Secondly, it is difficult to reconstruct the auto power 
spectrum of the target line f rom the cross spectrum of the 
two maps ilCong et al.l il2014l ~). hereafter G14). Thus we seek 
a method to obtain the auto spectrum of a single map with¬ 
out foreground contamination. 

Though “foreground” lines could technically come from 
lower or higher redshifts than the target line, only lower red- 
shift lines are likely to pose a problem, as the signal from an 
intensity mapping survey typically grows weaker with red- 
shift. The fact that these lines come from lower redshifts 
offers a potential way to remove them. At lower redshifts, 
galaxy masses tend to be larger, so we expect there to be 
more very bright sources of a foreground line than a tar¬ 
get line. This means that the brightest pixels in a survey 
will tend to be foreground galaxies, and the foreground con¬ 
tamination can be at least partially removed by masking 
out the br ightest pixels in a survey. Th is technique was dis¬ 
cussed by IVisbal. Trac. fc LoebI ll201ll ') , who found that it 
tends to bias the target power spectrum because some sig¬ 
nal is masked along with the foregrounds. However, as we 
will show below, much of the cosmological information in 
the power spectrum is preserved after masking, even though 
most of the astrophysical information is lost. G14 also ex¬ 
plored this technique for Lya, though they do not appear to 
show any biasing of the signal power spectrum. 

The aim of this work is to explore the effects of 
pixel masking on intensity mapping surveys in detail. 
We perform our study us ing simulated intensity maps. 
IVisbaL Trac. fc LoebI (l201ll ~) used N-body dark matter sim¬ 
ulations for their analysis, but this technique is more numer¬ 
ically intensive than required for our purposes. We therefore 
describe below a different method for simulating intensity 
maps, in which we assume a model for the matter power 
spectrum is known a priori. This makes it possible for us 
to quickly simulate relative large areas of the sky. We use 
the m atter power spectrum from GAME ilLewis fc Ghalliri^ 
[ 2 ^ along with some empirically estimated luminosity 
functions to generate signal and foreground intensity maps. 

Using these simulations, we will demonstrate below that 
pixel masking is an effective technique for removing fore¬ 
grounds when the signal is considerably brighter on average 
than the foregrounds, or when the pixel size of the survey 
is small enough that individual foreground sources can be 
isolated effectively. We find that GO intensity maps meet 
the first criterion, and Lya intensity maps presumably meet 
the second. GII maps on the other hand, likely will not meet 
either criterion, and therefore cannot be easily cleaned using 
simple pixel masking. 

It _ is _ important to note that 

iBrevsse. Kovetz. fc Kamionkow^ ||2014| ') (hereafter BKK) 


found that the amplitude of a GO intensity mapping power 
spectrum is extremely uncertain, and that different models 
yield very different results. We have every reason to believe 
that the modeling of each and every signal and foreground 
line we discuss here is similarly uncertain. Therefore it 
should be noted that the details of the luminosity function 
modeling described in this paper should be taken with a 
grain of salt and the exact amplitudes of the power spectra 
discussed below could vary significantly from the values 
we use. Nevertheless, our intention is to explore the effects 
of pixel masking on contaminated intensity maps, and our 
general conclusions should hold for many different models. 

This paper is organized as follows. Section 2 describes 
the method for calculating a power spectrum for an intensity 
map, then summarizes the models used for the three signal 
lines we consider. Our simulation methods are explained in 
Section 3, first for maps with no large scale structure, and 
then for maps with clustering included. Section 4 explains 
our method of removing foreground lines from these maps 
through pixel masking, and Section 5 shows the results of 
this masking. These results and pertaining issues are then 
discussed in Section 6, and we conclude in Section 7. 


2 SIGNAL & FOREGROUND MODELS 


For all of the lines we discuss in this paper it is reasonable 
to assume that all of the emission in a given line comes 
from within individual galaxies. We also assume that the 
galaxy luminosities in a given line are uncorrelated with one 
another. Under these assumptions, the three dimensional 
power spectrum of an intensity map takes a fairly simple 
form 

P{k, Z) = (Tf (2)Pgal(A:, Z) + Pshot(2). (2.1) 


The first term in this expression contains the sky-averaged 
brightness temperature (T) of the line, multiplied by the 
galaxy power spectrum Pgai, which is the linear dark mat¬ 
ter power spectrum multiplied by a bias factor (BKK). This 
term, which we refer to as the clustering term, gives the con¬ 
tribution to the power spectrum from the large scale struc¬ 
ture of matter in the universe. If the intensity field we mea¬ 
sure were spatially continuous this would be the entire power 
spectrum. However, since the signal comes from a large num¬ 
ber of randomly placed discrete sources, we must include an 
additional scale-independent shot noise term Pshot. 

The values of (T) and Pshot depend on the astrophysical 
conditions within the emitting galaxies. We quantify this 
astrophysical information through the luminosity function 
<E>(P), which gives the comoving number density of halos 
with luminosities in the desired line between L and L + dL. 
In terms of '1?(P), the average temperature and shot noise 
are 


(T) (z) = 

and 

-fshot(-2^) ~ 


1 (l + zf 


Stt 


'^em 


kbH{z', 

1 (1 + ^)^ 


/' 

Jl, 


Stt 


kbH{z) 


"f 


L^{L)dL, (2.2) 


P^$(P)dP, (2.3) 


where t'em is the emission frequency of the line, c is the 
speed of light, kb is Boltzmann’s constant, and H{z) is the 
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Hubble parameter. For some lines, the luminosity function 
can be difficult to measure, and it is often easier to estimate a 
relation L{M) between the luminosity of a halo and its mass. 
In this case, the integrals in equations (Ell) and (E3 are 
over mass instead of luminosity, with the integrands replaced 
by L{M)dn/dM and (M) dn/dM. In th is case, we use the 
mass function dn/dM from [Tinker et ^ (I 2 OO 8 I I instead of 
the luminosity function to estimate the number of galaxies 
with a given luminosity. 

Real intensity mapping surveys will survey the sky in 
several frequency bands, which corresponds to several differ¬ 
ent redshift slices. The intensities in these slices will be cor¬ 
related due to the existence of line-of-sight Fourier modes. 
However, we can obtain a reasonable approximation of a 
true survey by treating each frequency band as an inde¬ 
pendent map, which can then be stacked with the others to 
improve signal to noise. In order to facilitate the comparison 
of maps at different redshifts, we define our maps in angular 
coordinates and study the angular power spectrum Ci in a 
single slice. If we collapsed our whole 3D volume down to 
2D, then we would lose a large amount of the information in 
our map. Discussing the angular power spectrum of individ¬ 
ual slices, however, only sacrifices the information present in 
the line-of-si g ht m odes. This is the approximation used in 
IPuIjen et al.l (l2013l l and [Brevsse. Kovetz. fc Kamionkowsl^ 
1 2014) . All of the maps discussed in this paper should be 
thought of as a single slice of a full 3D survey. 

Since intensity mapping experiments are likely to have 
fairly narrow frequency bands, it is reasonable to assume 
that the quantity (T) does not change significantly over the 
width of a single band. This means that we can write the 
angular power spectrum as 

Ct{z) = {Tf {z)Cf^ + C'f°\ (2.4) 


The galaxy angular power spectrum is calculated from Pgai 
using 



k^P{k) 


f{r)je{kr)dr 


dk, 


(2.5) 


where r is the comoving distance, ji{kr) is the spherical 
Bessel function and /(r) is the selection function which is 
determined by the ?. For simplicity, we assume a top hat 
/(r). This integral is computationally difficult to evaluate, 
so we approxi mate it in the high-^ limit using the Limber 
approximation llLimbeill 19531 : lRubirJ[i954fl and in the low-^ 
limit by assuming that /(r) is a delta function. More details 
on these approximations can be found in Section 2 of BKK. 
The shot noise spectrum C|*'°‘ is scale independent, so the 
Limber approximation can be used on all scales. 

The values of (T) and Pshot for a given line depend 
sensitively on the exact shape of the luminosity function of 
that line. This luminosity function in turn depends on the 
conditions within the emitting galaxies which are highly un¬ 
certain, especially for high redshift sources. We can attempt 
to predict L{M) or 'l'(L) for different lines using various 
empirical observations. However, BKK found that differ¬ 
ent models assumed for CO emission yielded power spectra 
which spanned roughly two orders of magnitude in ampli¬ 
tude. Therefore the models we summarize below should not 
be interpreted as precise predictions of the power spectra for 
these lines. Rather, they are intended as a means to gain un¬ 
derstanding of how the shapes of the luminosity functions of 


foreground lines can affect how strongly they contaminate a 
target line, separately from the overall amplitude ambiguity. 


2.1 CO 


When modeling the emission of CO and its foreground lines, 
we assume that only halos with masses above some cutoff 
mass Mmin (assumed here to be 10® Mq) can emit the line in 
question, and we assume that only a fraction /duty equal to 
the timescale of star formation over the age of the universe of 
halos are emitting the line at any given time. For generality, 
we assume that the luminosity of a halo is a power law in 
its mass 

= A/duty , (2.6) 

where A and b are f r ee par ameters. _ _ 

In iPullen et al.l ll2013ll and iLidz et ^ l^2011^ , the cal¬ 
culation of L{M) for CO uses a series of empirical scaling 
relationships, starting with a relation between the fine lumi¬ 
nosity and FIR luminosity of the form 


r / r' \ 

TfIR _ „ j Tiine \ 

Lq “ Vk km s-1 pc^J 


(2.7) 


where Cfir and Xfir are constants set through observa¬ 
tions. In the above relation, the line luminosity is given in 
units commonly used for spectral line observations. The con¬ 
version to solar luminosities is given by 


-^line _ o y \ ^ _ -^line _ 

Lq^ U GHz/ Kkms-ipc2’ 


( 2 . 8 ) 


llCarillill201lll . These relations, in combination with the star 
formation rate-FIR lumino sity and halo mass -star forma¬ 
tion rate relations given in iPullen et al.l ll2013fl . we get the 
following expressions for A and b: 


M = 3 X 10”“ 


/6.5X 

V Cfir ) V1 GHz 


(2.9) 


b = (2.10) 

For GO, we follow BKK and IPullen et al.l (I 2 OI 3 I I in using the 
IWang et al.l 1 I 2 OIOII GO-FIR relation with Cfir = 1.4 x 10“® 
and -Yfir = 5/3. This g ives t he values A = 2 x 10~® and 
6=1 from IPullen et aP (I 2 OI 3 II . 

When considering possible foreground lines, we consider 
only lines with lower rest frame frequencies than our target 
lines. This is because, as shown in G14, projecting two power 
spectra from different redshifts boosts the amplitude of the 
lower redshift spectrum relative to the higher redshift one. 
Because we are working here with angular power spectra 
instead of three dimensional power spectra, this projection 
effect is included naturally in our calculations. For a GO sur¬ 
vey targeted at 2 = 3, we are therefore concerned with lines 
that have emission frequencies between the 115 GHz CO 
rest frame frequency and the 28.8 GHz observing frequency. 

The CO (1-0) line is expected to be considerably brighter 
on average than any other line in this range, i.e. (Too) 
should be much greater than (T) for any foreground line. 
If we assumed that L{M) were linear in M for all lines this 
would mean that no foreground line could dominate over the 
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CO signal. However, the addition of the power law depen¬ 
dence on M means that the CO map could still be contam¬ 
inated by the brightest sources for a foreground line which 
has fe > 1. The shot noise term in particular is very sensitive 
to the value of b. For more general luminosity functions than 
the one we use here, this would mean that a line is poten¬ 
tially a problem if it falls off more slowly at high luminosities 
than the target line. 

We consider foreground lines emitted by four molecules: 
HCN, HCO+, CN, and CS. All of these molecules have 
higher critical densities than CO, and thus tend to trace 
denser regions of galaxies. For each molecule, we use an em¬ 
pirical correlation with FIR to estimate A and b. The results 
and sources for each line are given in Table [T] along with 
those of CO for reference. Note that for the CN and CS lines 
the FIR relations used were measured for higher order tran¬ 
sitions (3-2 and 7-6, respectively). Thus we have made the 
assumption that the intensities of these lines are indepen¬ 
dent of which transition is being considered. The situation 
in reality is likely not so simple, but this approximation will 
suffice for our purposes. 

Figure[T] shows the clustering and shot noise power spec¬ 
tra for each of these lines compared with that of CO. Because 
all of the foreground lines have b > 6co, they are all domi¬ 
nated by shot noise. A possible qualitative reason for this is 
the fact that all of these lines tend to trace denser gas than 
CO, and thus their luminosities may be more sensitive to 
the environments in their host galaxies. 

Of the lines we consider, all but HCN have power spec¬ 
tra well below that of CO. The HCN line however, actually 
starts to dominate over CO at ^ ~ 1000. Though HCN is 
only ~ 4% as bright as CO on average, it has a higher value 
of h, which means that it produces a small number of bright 
sources which contribute a large amount of shot noise. This 
demonstrates how a line with relatively small average inten¬ 
sity can still be a problematic foreground. For the remainder 
of this paper, we will consider only the HCN foreground line 
since in our modeling the other lines are subdominant. 

The shaded region in Figure [T] shows a very crude esti¬ 
mate of the theoretical uncertainty in the HCN power spec¬ 
trum, roughly an order of magnitude. BKK found that the 
amplitude of the CO spectrum could also vary by roughly an 
order of magnitude in either direction. For convenience we 
only show this shading for HCN, but all of the other spectra 
are similarly uncertain, if not more. Given the vast amount 
of uncertainty both in the modeling of these spectra as well 
as in the empirical measurements used in our model, we do 
not attempt to make a more accurate estimate of the error 
bars on these spectra. 


2.2 Lya 

For the Lyman a line we consider a hypothetical survey 
targeted at 2: = 7, and we follow the modeling of G14 and 
references therein, which we summarize briefly here. Since 
the Lya line will be observed in the infrared rather than 
the radio, we use intensity units here instead of brightness 
temperature. Recombination and collision processes within 



Figure 1. Power spectra for lines described in Table [T] The solid 
blue line is the target CO power spectrum, and the dashed lines 
show the spectra for the foreground lines. Note that all of the 
foregrounds are shot noise dominated, and only HCN is compa¬ 
rable to CO. The shaded region shows a rough estimate of the 
theoretical uncertainty in these power spectra. 


galaxies give a Lya luminosity 


= 5.1 X 10®/Lyc(2) [1 - /‘“(M)] 


SFR{M, z) 
Mq yr-i ’ 
( 2 . 11 ) 


where fhya = 3.34 x 10“® x (1 -|- 2)^'®^ is the frac¬ 
tion of Lya photons not absorbed by dust, and /e™ = 
exp (—5.18 X 10“® X is the escape fraction of ion¬ 

izing photons. The star formation rate is parameterized by 


SFR{M,z 


7) 


X 


1.6x10” 


M 


1-b 


M 

Ah 


1-b 


2.59 

1 

^ Ml 

Ad ® 

-2.25 

\ 

Ads j 

1 


^ -0.62 
( 2 . 12 ) 


We take Admin ~ 10® AIq and Mmax = AIq. Note that 
for simplicity we neglect Lya emission from the IGM. The 
IGM contribution is small compared to the halo emission, 
and it is not as easily simulated using our methods. We 
calculate the bias for Lya following G14. 

For Lya at 2 = 7, three foreground lines are considered: 
Ha coming from 2 ~ 0.5, OIII from 2 ~ 0.9, and Oil from 
2 ~ 1.6. For these lines, the luminosity function is assumed 
to be a Schechter function 

$(L)dL = exp(^-^^^, (2.13) 


where 'I>(L) is the comoving number density of halos with 
luminosities between L and L + dL, and St, L,, and a are 
parameters which we obtain from iLv et al.l (l2007l L With the 
above relations, we can calculate angular power spectra for 
Lya and the three foreground lines using equation (12:4)1 . 
Following G14, we assume that the halo bias for these lines 
is proportional to the halo mass. Figure 0 shows the calcu¬ 
lated power spectra for Lya and the three foreground lines 
along with the total foreground spectrum. Note that the 
foreground lines dominate entirely over the signal, and un- 


© 0000 RAS, MNRAS 000, 000-000 


















Masking line foregrounds in intensity mapping surveys 5 


Table 1. Various parameters for the target CO(l-O) line as well as several possible foreground lines. Parameters include the emission 
frequency fem, the parameters of the FIR correlation CpiR and Vpijr, the L{M) parameters A and b, and the observable parameters 
(T) and 


Line 

k'em (GHz) 

Cfir 

IfpiR 

A 

b 

(T) (/rK) 

Cf°* 

Source 

CO(l-O) 

115 

1.35x10“® 

1.67 

2x10“® 

1 

0.60 

7.8x10“’^ 

P13 

HCN(l-O) 

88 

794 

1.0 

1.7xl0“i® 

1.67 

0.023 

2.1x10“® 

Gao & Solomon 120041 

CN(l-O) 

113 

1.6x10'^ 

0.89 

6.9xl0“i® 

1.87 

0.010 

1.0x10“'^ 

Riechers et al. f2007') 

CS(l-O) 

49 

2.1x10"^ 

1.0 

2.7xl0“i® 

1.67 

2.1x10“"^ 

3.8x10“® 

Zhane et al. ('2014') 

HCO+ 

89 

158 

1.11 

1.2xl0“i® 

1.5 

1.94xl0“'‘ 

3.4x10““ 

Riechers et al. ('2006') 




Figure 2. Theoretical power spectra for Lya and its three fore- Figure 3. Theoretical power spectra for CII and its four fore¬ 
grounds as well as the total foreground spectrum ground CO lines along with the total foreground spectrum. 


like CO, the foregrounds have a signihcant clustering com¬ 
ponent. 


2.3 CII 


We base our simulations of CII intensity maps on the mod¬ 
eling done by I Silva et al.l ll20l4i . hereafter S14. We consider 
a possible survey targeted at z = 7. For the CII luminosity 
we use model m2 from Table 1 of S14, where the luminosity 
is given by 


log 10 


u© 


log 10 


/ tp{M,z) \ 

\ Mq J 


+ 6.9647, 


where 


ij{M,z) = Mo 


\Ma) 



(2.14) 


(2.15) 


The parameters in the formula for ip at z = 7 are Mq = 
6.6 X 10"® Mq, Ma = 10® Mq, Mt, = 1.6 x 10“ Mq, 
ttcii = 2.25, and bcii = —2.3. The average CII intensity and 
shot noise can then be calculated using equations (EH) and 
(EH above. 

The primary foregrounds for a CII survey come from 
CO. Since such a survey would observe at roughly 240 GHz, 
we need not worry about the CO (1-0) line at 115 GHz or the 
CO(2-0) line at 230 GHz. The higher order transitions are 
potentially problematic though. We estimate the intensities 
of the higher order CO lines using our formulae above along 


with the line ratios given for su bmillimeter galaxies in Ta¬ 
ble 2 of ICarilli fc Walt^ ll2013li . This table gives the ratio 
^V^co(i-o) ^he transitions from CO(2-l) to CO(5-4). 
Following S14, we consider the CO(6-5) line as well, and as¬ 
sume that it has the same luminosity ratio as the CO (5-4) 
transition. We also calculate the bias by assuming that the 
halo bias is proportional to the halo mass. Figure |3] shows 
the calculated power spectra for these CO lines as well as the 
spectrum of the target CII line. As with Lya, the foreground 
lines are strongly clustered and considerably brighter than 
the target line. 


3 SIMULATIONS 

In order to explore the effects of foreground contamination 
in more detail it is useful to create simulated maps of var¬ 
ious lines. In this section we describe our simulation meth¬ 
ods, first for maps with shot noise power spectra, then for 
maps with large scale clustering included. We will explain 
our simulation process using CO at z = 3 as an example. 


3.1 Shot noise maps 

When simulating a map, we first define a grid of pixels, 
with the solid angle of each pixel defined by the angular 
resolution of a hypothetical instrument, the pixel depth by 
the spectral resolution, and the total area dehned by the 
proposed survey area. For our CO simulations, we use a 
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fiducial survey observing at 30 GHz (corresponding to CO 
at 2 = 3) covering 550 sq. deg. with an angular resolution of 
10 arcmin and dv /v = 10~^. These n umbers are similar to 
those discussed in Pullen et al. and BKK. Since each 

pixel in this simulated map is large compared to a galaxy, 
we can assume that the line emitters are essentially point 
sources. The expected mean number of sources in each pixel 
is given by 

{N) iz) = nPpix/duty = Hpix(«)/duty(«) H ^^dM, 

-'Mco.min 

(3.1) 

where n is the total halo number density and Ppix is the 
comoving volume of a single pixel. 

With no clustering taken into account, galaxies should 
be Poisson distributed on the sky. Thus we can draw the 
number of sources in each pixel from a Poisson distribution 
with mean (N). Assuming the line luminosities of individual 
sources are uncorrelated, we can then randomly assign a 
luminosity to each source. Since we are using a simple model 
in which luminosity is determined by mass, we first draw a 
mass from the distribution 

then calculate the fine luminosity from L{M) for each galaxy 
and convert this to a brightness temperature to get the fin¬ 
ished map. 

Since we are simulating relatively small regions of the 
sky, we can calculate power spectra in t he fiat sky ap¬ 
proxim ation. similar to the method used in IChiang fc CherJ 
(l2012ll . For simplicity, we assume the survey area is a square 
located near the equator of whatever sky coordinates are in 
use. 

Consider a map with Apix pixels at positions x = 
{xi,yi), where Xi and j/i run from 1 to ^Npi^- The intensity 
at each pixel r(x) can be decomposed into Fourier modes 
Ok through the discrete Fourier transform, 

Uk = ^ V . (3.3) 

-/Vpix 

X 

The angular power spectrum of the map is then 

^£=27rfc/v/n ~ HCfc, (3.4) 

where Q is the total solid angle of the survey and 

Ck = {W\^). (3.5) 

The average in equation (I3.5|l is taken over all k modes in the 
interval k—1/2 ^ |k| < A:-1-1/2. Note that in this calculation 
X and k are in units of pixels and pixels”^ respectively. 

The left hand column of Figure |4] shows a map of CO 
emission simulated using the above technique as well as 
its power spectrum. The power spectrum has the expected 
scale-independent form of a shot-noise dominated sample. 
The solid line is the predicted shot noise amplitude calcu¬ 
lated using equation (lOl) for CO at 2 = 3, included to show 
the consistency between our model and our simulations. The 
scatter of the points around this fine is cosmic variance error 
due to the finite size of our simulation. There is no instru¬ 
mental noise included in the simulated map. 


3.2 Adding clustering 

The simulations above assumed that galaxies are randomly 
distributed on the sky, but in reality the large scale struc¬ 
ture of the universe imposes a pattern of clustering on 
the galaxy distribution. The typical method of simulating 
maps with realistic clustering is to perform an N-body dark 
matter simulation and pop ulate the result with galaxies 
llVisbal. Trac. fc LoebluOlUl . However, this is computation¬ 
ally intensive and makes it difficult to simulate large sky 
areas. Since we have predictions for the angular power spec¬ 
tra of these maps we choose to generate random density 
fields with the desired power spectra rather than attempt to 
simulate the entire history of the region being simulated. 

We want to simulate a density field (5(x) = (n(x) — h)/h 
which includes clustering. The 5 map should have zero mean 
and its power spectrum should be the halo power spectrum. 
In a pixel located at x, instead of drawing the number of 
galaxies from a Poisson distribution with mean {N ), we use 
{N) {5{:k) + 1). Since galaxy luminosities do not depend on 
the clustering, the procedure for assigning luminosities is the 
same as before. This will give a map with a power spectrum 
that contains both clustering and shot noise components. 

It only remains then to generate a density field with the 
desired halo power spectrum. If the field is Gaussian, the 
process is fairly straightforward. In a full-sky Gaussian ran¬ 
dom field, each spherical harmonic coefficient airn is drawn 
from a Gaussian distribution with variance %/C/. In the flat 
sky approximation, we draw each Fourier mode with magni¬ 
tude in the range fc — 1/2^ |k| <fc-|-l/2 randomly from a 
Gaussian with variance Ck- This map can then be converted 
back to spatial coordinates using the inverse discrete Fourier 
transform. One important thing to note is that the resulting 
map must be real, while the Fourier modes are complex. This 
means that we must impose the condition that a_k = a^- 

Unfortunately, at the redshifts we are considering, the 
galaxy distribution is highly non-Gaussian. Attempting to 
impose a Gaussian distributed S(x) produces pixels where 
(N) (S(x) + 1) is negative, which is obviously unphysical. A 
better approximation for the ga laxy distribution wo uld be 
to use a log-normal distribution (IColes &: Jonesl[l99lll . This 
distribution has the important property that it is zero for 
negative densities and it appears to be a reasonably good 
fit to the observed galaxy distribution. The procedure for 
generating a log-normal random field uses the fact that that 
a log-normal map 5ln{x.) can be generated from a Gaussian 
map with variance ctq using 

5ln{x) = ^ (3 g ) 

There exists a convenient relation between the correlation 
functions ^(r) of two maps related in this way 

5G(r) = In [1-bCniv(r)] . (3.7) 

This allows us to generate a log-normal random field Sln 
with the desired characteristics. 

We start with the galaxy power spectrum and con¬ 
vert it to a flat sky approximation From this, we can 

calculate the correlation function we want our log-normal 
field to have using 
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Figure 4. (Top left) Simulated shot noise dominated CO intensity map from 2 = 3. (Bottom left) Power spectrum of simulated map 
(red points) overplotted with the predicted shot noise power spectrum from equation II2.3II . (Top right) Simulated CO intensity map 
from 2 = 3 with clustering included. (Bottom right) Power spectrum of simulated clustering map overplotted with the predicted power 
spectrum including both clustering and shot noise components from equation 


We then calculate using equation (ITfl) and convert this 
to a new angular power spectrum using 

2 

Ck ~T r^G(r) sin(fcr)dr. (3.9) 

k Jo 

We can then draw a Gaussian random field So with this 
power spectrum using the procedure outlined above, then 
convert it to a log-normal random field using equation (EH). 

The result of this process is a map like the one shown 
in the top right panel of Figured with a power spectrum as 
shown in the bottom right panel. It is obvious when com¬ 
paring these plots to those in the left hand column that the 
second map is much more strongly clustered and that the 
power spectrum contains a distinct scale-dependent compo¬ 
nent. The blue curve in this plot shows the predicted CO 
power spectrum from Section 2.1. This predicted spectrum 
is in good agreement with the simulated data, which again 
shows the consistency between our theoretical models and 
our simulations. 

3.3 Simulation results 

With this method in place, we can now simulate maps for 
any line for which we can define either a function L{M) 
or a luminosity function $(i). For lines where we use a 
luminosity function, we simply replace the mass function in 
equation ES with ^(L). Since we are working in angular 
space, it is also very easy to combine maps for different lines. 


A map containing two lines is simply the sum of the maps 
for the individual lines. If we simulate a map for CO at 
2 = 3 and a foreground HCN line at 2 = 2, we get power 
spectra like the ones shown in Figure [5] The solid curves in 
this plot are the best fit power spectra obtained by treating 
(T) and Cf as free parameters. As expected, the effect of 
the foreground contamination is to significantly increase the 
shot noise of the map, which is particularly obvious at high 
£. 

The amplitudes of the power spectra in these simula¬ 
tions depend on exactly how many halos are drawn and ex¬ 
actly what masses are assigned to them. This means that the 
result can vary somewhat from the theoretical predictions. 
The HCN power spectrum is particularly sensitive to this, 
since it depends so much on the highest mass halos. The 
shaded region in Figure El shows the 95% confidence range 
of the HCN power spectrum, calculated by comparing the 
results of 500 simulated HCN maps. 

Figure El shows simulated maps for CO at 2 = 3 and 
Lya and CH at 2 = 7 along with their foregrounds in order 
to illustrate the qualitative differences between them. The 
CO simulations cover 550 deg^ with 10 arcmin resolution 
and Aiz/i'obs = 10“^, the Lya simulations cover 1 deg^ with 
0.1 arcmin resolution and Aiz/i/obs = 1/40 (G14), and the 
CH simulations cover 100 deg^ with 3.2 arcmin resolution 
and Av/vahs = 1.7 x 10“® (S14). 

The HCN foreground map is considerably fainter on 
average than the CO map, but though the brightest sources 
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Figure 6. Simulated maps of CO at 2 = 3, Lyo at 2 = 7, and CII at 2 = 7 (top row, left to right) along with their foregrounds (bottom 
row). The CO simulations cover 550 deg^ with 10 arcmin resolution, the Lya simulations cover 1 deg^ with 0.1 arcmin resolution, and 
the CII simulations cover 100 deg^ with 3.2 arcmin resolution. 


are difficult to see on this image the colorbar shows that the 
HCN map extends to higher intensities than the CO map. 
The Lya maps cover a much smaller region of the sky than 
the others, so the clustering features in this image appear 
much larger. In addition, each pixel in the Lya foreground 
map covers a rather small volume of space, so most of the 
pixels in this map are dark. Though the Lya foregrounds are 
brighter than Lya on average, most of the intensity in this 
map comes from a small number of pixels. The foreground 
map for CII is much brighter than the signal, and the pixels 
in the CII map are large enough that there is strong emission 
in most pixels. 

Though the models we use in this work are fairly sim¬ 
ple, this simulation method is straightforward to generalize 
to different, more complex models. Any power spectrum and 
luminosity function can be input to get a simulated map. 
One important effect that we have not taken into account 
here is the correlation between maps made at different fre¬ 
quencies due to line-of-sight Fourier modes. Since the target 
lines we consider here are widely separated from the fore¬ 
ground fines in frequency space this effect should not be 
significant for the problems discussed here. However, it will 
need to be taken into account when attempting to accurately 
simulate the three-dimensional data taken by realistic inten¬ 
sity mapping experiments. 


4 METHOD 

As noted above, the majority of the contamination in CO 
surveys comes from a few bright foreground emitters which 
add a large amount of shot noise to a survey. As noted by 
G14, this means that the foreground effect could be mit¬ 
igated if we simply mask out the brightest pixels in our 
survey. Figure [7] shows the number of sources/pixel which 
produce brightness temperatures above a given value. At 
brightness temperatures above a few hundred K, there are 
more HCN emitters than CO emitters. This is primarily due 
to the fact that 6hcn > 6co, with some additional contri¬ 
bution from the fact that there are more high mass halos at 
2 ~ 2 than at 2 ~ 3. 

This means that if a given pixel has a very high inten¬ 
sity, it is likely that the extra flux is coming from a fore¬ 
ground HCN emitter. Thus if we mask all pixels brighter 
than a given value we will remove on average much more 
foreground emitters than target emitters. In addition, the 
foreground sources we mask are exactly the bright sources 
which produce the worst contamination. Therefore we ex¬ 
pect that we can clean the foregrounds out of a map by 
masking out the brightest pixels. 

The effects of pixel masking on the power spectra of 
these maps are difficult to understand analytically. We can 
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Figure 5. Power spectra of simulated maps of CO at z = 3 
(blue), HCN at 2 = 2 (green), and the sum of the two (red). Solid 
curves are best fit power spectra allowing (T) and Cf to vary. As 
expected, the HCN foreground contributes a significant amount 
of extra shot noise at high £. The shaded region corresponds to 
the 95% confidence interval for the simulated HCN maps. Varia¬ 
tions from simulation to simulation can cause the power spectrum 
amplitude to change significantly. 



Figure 7. Number of sources per pixel which contribute bright¬ 
ness temperature greater than T. At low intensities, the target 
CO line dominates as expected. However, at high intensities the 
HCN emitters begin to dominate because the HCN luminosity 
function does not fall off as quickly as that of CO. 

make a very rough estimate if we assume that we can per¬ 
fectly isolate and remove all sources which contribute an 
intensity greater than some cutoff. The ratio of the average 
brightness temperature of a line before and after applying 
this cutoff is 

unmasked 

If we choose to remove ~ 1% of the pixels from a map, Mcut 
is the mass of a halo such that 1% of halos have M > Mcut- 
For CO, we have Mcut = 6.2 x 10^° Mq and for HCN we have 
9.5 X 10^° ^Q- The ratio of the shot noise before and after 


masking is similar to equation (HU with L{M) replaced 
with L^{M). 

This analytical procedure predicts that the ratio of the 
masked to unmasked CO intensity will be ~ 0.7 and that of 
the masked to unmasked CO shot noise will be ~ 0.06. The 
ratio of masked to unmasked HCN shot noise is predicted 
to be ~ 2 X 10“®. Therefore, although masking pixels causes 
the CO power spectrum to decrease, we expect the HCN 
spectrum to decrease far more, leaving behind a signal dom¬ 
inated power spectrum. However, the procedure described 
here neglects a few key effects which would be present in a 
real map. If a map pixel contains a bright source, masking 
it will also remove all of the fainter sources present in the 
same pixel. In addition, some pixels will be masked because 
they contain a large number of faint sources rather than a 
single bright one. In order to fully account for these effects, 
we need to rely on our simulated maps, as we do below. 

One question that arises naturally is exactly how many 
pixels should be masked? For CO, and any other hypothet¬ 
ical line where all of the foregrounds are shot-noise domi¬ 
nated, the answer is simple. We only need to mask until the 
masked power spectrum shows clustering behavior out to the 
desired angular scale. Even if the CO foregrounds are faint 
enough that CO dominates entirely, this masking will reveal 
the CO clustering behavior on scales which are normally ob¬ 
scured by shot noise. For lines with clustered foregrounds, 
the answer is less clear. If we have a reasonable estimate of 
the luminosity functions of the signal and foregrounds, we 
can use a plot like Figure[7]to predict a cutoff intensity (this 
is the method used in G14). 


5 RESULTS 

Here we show the results of applying the masking procedure 
described above to CO, Lya, and CH intensity maps. 

5.1 CO 

Figure [8] shows the effect of pixel masking on the power 
spectra of three simulated maps: one with just CO, one with 
HCN, and the sum of the two. In order to make the effects 
of masking more obvious, we have simulated the HCN map 
with the value of the A parameter increased by a factor of 3. 
This creates a map where the foreground shot noise entirely 
dominates over the signal. Given the uncertainties in the 
power spectrum modeling, it is not impossible for this to be 
the case in reality. 

In Figure[8l the dots and dashed curves show the power 
spectra and best fit curves for the three maps with no mask¬ 
ing. The boosted foreground power spectrum entirely dom¬ 
inates on all scales. The pluses and solid curves show the 
spectra of the maps after all pixels brighter than 7 fiK are 
masked. This value corresponds to masking roughly 1% of 
the pixels in the map with both signal and foregrounds. Af¬ 
ter masking, the foreground power spectrum has dropped 
dramatically, and the CO power spectrum has fallen by a 
much smaller amount. But most importantly, the red total 
power spectrum is very similar to the CO power spectrum. 
The power spectrum of the map with foregrounds is nearly 
identical to that of a map without them. There remains some 
small amount of shot noise contamination, but this could be 
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Figure 8. Power spectra and best fit curves for simulated maps 
of CO at redshift 3 (blue), HCN at redshift 2 (green) with the 
amplitude boosted to dominate over CO, and the sum of the two 
(red). Dashed curves/dots show the spectra before masking, solid 
curves/pluses show the results of masking all pixels above 7 pK 
(~ 1% of the pixels in the total map). After masking, the maps 
with and without foregrounds have very similar power spectra. 

removed by choosing a lower cutoff value for masking. Thus 
it appears that foreground contamination in CO intensity 
maps can indeed be mitigated by masking bright pixels. 

It is worth mentioning that the amount by which the 
CO power spectrum drops differs somewhat from the sim¬ 
ple prediction made in Section 4. The clustering term in our 
simulation decreases roughly 30% less than the simple cal¬ 
culation predicts, and the drop in shot noise is ten times less 
than predicted. As mentioned above, this happens because 
the calculation in Section 4 does not take into account the 
fact that there are multiple sources in each pixel, which will 
either be masked along with a single bright source or add 
together to mimic a bright source. This discrepancy makes 
it difficult to use equation (4.1) to predict the unmasked 
spectrum from a masked one. 

5.2 Lya 

To facilitate comparison between our results and those given 
in G14, we choose to mask 3% of the pixels in our Lya sim¬ 
ulations. Figure [9] shows the effectiveness of this masking. 
The dashed lines show the power spectra of the signal, fore¬ 
grounds, and total map before masking, and the solid lines 
show the spectra after masking. For the sake of visibility, we 
have plotted only the total foreground spectra rather than 
those of the three individual lines. After masking, the fore¬ 
ground has dropped below the signal and the map which 
includes both signal and foregrounds gives a spectrum very 
similar to the map without foregrounds. Thus it appears 
that for the model and telescope resolution simulated here 
it is possible to remove most of the foreground contribution 
to the power spectrum. 

5.3 CII 

Figure [To] shows the power spectra of the CII simulations be¬ 
fore and after masking 1% of the pixels. Note that unlike in 



Figure 9. Power spectra and best fit curves for simulated maps of 
Lya at 2 = 7 (blue), the sum of the three foreground lines (green), 
and all four lines combined (red). Dashed curves/dots show the 
spectra before masking, solid curves/pluses show the results of 
masking all pixels above ~ 100 Jy/sr (3% of the pixels in the 
total map). As with CO, the masking removes a large amount of 
the foreground contamination. 



Figure 10. Power spectra and best fit curves for simulated maps 
of CII at 2 = 7 (blue), the total CO foreground (green), and 
the total emission (red). Dashed curves/dots show the spectra 
before masking, solid curves/pluses show the results after masking 
all pixels brighter than ~ 10^ Jy/sr (1% of the pixels in the 
total map). After masking, the foregrounds still dominate over 
the signal. The foregrounds continue to dominate no matter what 
masking percentage is chosen. 


our CO and Lya simulations, the foreground lines dominate 
over the signal even after masking. Though we have plotted 
the results of masking 1% of the pixels, the same basic re¬ 
sult holds true no matter how many pixels we mask. At this 
resolution, it is not possible to move the foreground power 
spectrum below that of the signal no matter how many pix¬ 
els are removed. 
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6 DISCUSSION 

The details of masking the simulated maps explain why the 
three lines we consider behave somewhat differently under 
masking. For the CO map, the contamination is shot noise 
from the few brightest galaxies, which is easily removed. 
However, both Lya and CII have foregrounds which include 
strong clustering components. In addition, for both of these 
lines the average intensity of the foregrounds is large com¬ 
pared to that of the signal, as seen in G14 and S14. Despite 
this, the masking works well for Lya and poorly for CII. 
The reason for this has to do with the size of the pixels used 
in each map. Lya is a higher frequency line than CII or 
CO, so it is easier to map with high angular resolution. This 
means that it is easier to isolate a single bright foreground 
source without removing as much signal. Our results show 
that if a CII survey is limited to arcminute scale pixels then 
foregrounds cannot be easily cleaned through pixel masking, 
and we must resort to cross-correlations to isolate the target 
line. 

The drop in the amplitudes of the target spectra af¬ 
ter masking is an unfortunate side effect of pixel masking. 
Ideally, an intensity mapping survey would be able to re¬ 
cover the dashed blue curve in Figure [8] rather than the 
solid blue curve. This is because we cannot avoid mask- 
ing some CO source s whe n removing HCN. As noted in 
IVisbal. Trac. fc LoebI (l201lll . some of the information in the 
power spectrum is lost in this masking process, specifically 
that encapsulated in the amplitudes of the clustering and 
shot noise terms. On a brighter note, since we are only mask¬ 
ing a percent or so of the pixels, the shape of the clustering 
term of the power spectrum does not change significantly 
when masked. This means that masking allows us to recover 
the shape of the galaxy power spectrum even on scales where 
the foregrounds dominate over the signal. 

This is a potentially useful meas urement for a variety of 
cosmo logical purposes. For example, iBrax. Clesse. fc Davis! 
(l2013tl show that modified gravity models alter the shape of 
the power spectrum which would be measured by a 21 cm 
survey. Similar changes could be studied in the masked maps 
we describe here. However, all of the astrophysical informa¬ 
tion which was contained in the amplitudes of the clustering 
and shot noise components is lost in masking, so we cannot 
make statements about the luminosity functions of the tar¬ 
get galaxies. This means that, while these intensity mapping 
surveys could still be useful with only this simple foreground 
cleaning, it may be necessary to use other foreground clean¬ 
ing methods to reach their full potential. 

Nonetheless, the information lost in pixel masking could 
potentially be recovered in a number of ways. Unfortunately, 
the estimate of the change in power given in equation (4.1) 
differs significantly from our simulations, so this cannot be 
used to recover the unmasked power. However, cross cor¬ 
relations between_different_jnaps_axe always possible, and 
S14 and IVisbal. Trac. &i LoebI ll201lll discuss the possibil¬ 
ity of using information in another wavelength to isolate 
pixels which contain foreground galaxies, thus limiting the 
amount of signal lost in masking. G14 also state that the 
behavior of line-of-sight Fourier modes may be different 
for signal and foreground lines. Another way could be to 
fit the one-point statistics of the m ap using P(D) analysis 
jLee. Ando, fc KamionkowskillioO^ l as the bright pixels are 


progressively removed, and compare the results with theo¬ 
retical simulations to retrieve the original amplitudes. This 
could be a means of obtaining the astrophysical information 
lost in masking (Breysse, Kovetz, & Kamionkowski 2015). 
Both the power anisotropy mentioned in G14 and the one 
point statistics could also be used to determine if the signal 
or a low-redshift foreground line dominates after a map has 
been masked. 

Note that care must be taken when fitting uncleaned in¬ 
tensity mapping power spectra to determine cosmological or 
astrophysical parameters. For example. The power spectrum 
of a GO intensity map can in principle be fit to determine 
the values of various model parameters (Li et al. 2015), and 
the fitted parameters will depend on how much shot noise 
is present in the map. However, the foreground line will add 
an uncertain amount of extra shot noise to the spectrum, 
which will bias the results of any fit. 

We have not included the effects of instrumental noise in 
any of these calculations. Though we leave for future work a 
full analysis of the behavior of noisy maps with masked pix¬ 
els, we can make some basic arguments to predict whether 
or not our basic premise would hold in a map with noise. If 
a map has too much noise, then the pixels which are masked 
will be bright due to random noise fluctuations in no corre¬ 
lation with the brightness of the target or foreground galax¬ 
ies in those pixels. Therefore, if a simulation has too much 
noise this pixel masking technique is useless. The survey 
parameters used in BKK, which were chosen to provide a 
reasonable chance of detecting CO at z = 3, give a noise 
per pixel aN = 1.7 pK. If we generate noise maps with the 
value in each pixel drawn from a Gaussian with zero mean 
and standard deviation of crjv, we can estimate how many of 
the masked pixels will be bright due to noise rather than sig¬ 
nal. We find that adding noise to our simulated CO-I-HGN 
map increases the number of pixels above our 7 pK cutoff 
by roughly 25%. These spurious bright pixels will reduce the 
effectiveness of the masking, but it should still be possible to 
remove most of the bright sources from a map. If the noise is 
significantly stronger, the masking will be ineffective. How¬ 
ever, in this case the signal to noise ratio for detecting the 
GO line at all becomes considerably smaller as well. 

In addition, we assumed the linear form for the under¬ 
lying matter power spectrum in our calculations. Though we 
account for some nonlinearity by using a lognormal galaxy 
density field, a full treatment of the nonlinearities would 
add more power to our maps on small scales. However, this 
would likely not have a significant qualitative effect on our 
results. There may be some minor differences in the ratio of 
signal/foreground power after masking, however the scales 
where nonlinearity becomes most significant are also scales 
where the power spectra tend to be shot noise dominated, 
so we do not expect any dramatic effects from a nonlinear 
calculation. 


7 CONCLUSION 

We have presented an exploration of the effectiveness of 
bright pixel masking on removing foreground lines from 
intensity maps. Using empirical luminosity function mod¬ 
els and simulated intensity maps we have illustrated how 
masking changes the power spectra of maps for three cases: 
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CO contaminated with HCN, Lya contaminated with var- 
ions atomic lines, and CII contaminated with higher order 
CO lines. For the CO survey, the foreground line was faint 
enough on average that removing the brightest pixels sig¬ 
nificantly dropped the amplitude of the foreground spec¬ 
trum. The high angular resolution possible in the Lya sur¬ 
vey meant that the foreground contamination was limited to 
a few pixels which could be easily masked. The CII survey, 
however, had both bright foregrounds and large pixels, so 
the masking was found to be ineffective. 

For all of the lines, masking bright pixels altered the am¬ 
plitude of the recovered power spectrum away from the de¬ 
sired uncontaminated value. This means that masking loses 
some of the information in the spectrum. However, in the 
two surveys where masking was effective, the masked spec¬ 
trum had a clustering component with the same shape as 
the unmasked clustering spectrum. Therefore, though the 
astrophysical content of the map is lost, the cosmological in¬ 
formation contained in the shape of the clustering spectrum 
can be recovered from a masked map. Thus, pixel mask¬ 
ing seems to be a useful technique for obtaining information 
from even a highly contaminated CO or Lya map. If we are 
to obtain the remainder of the information in these surveys, 
it will be necessary to use some other foreground cleaning 
technique, such as cross correlation, or to augment it with 
a P{D) analysis of the progressively masked power spectra. 
If we are to fully unlock all of the benefits of intensity map¬ 
ping surveys, it is imperative that we utilize these or other 
methods to isolate the signal from the foregrounds. 
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